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^_i I Abstract 

Higher order nets were introduced by Dick as a generalisation of classical (t, m, s)- 

nets, which are point sets frequently used in quasi-Monte Carlo integration algo- 
^^ ' rithms. Essential tools in finding such point sets of high quality are propagation 

rules, which make it possible to generate new higher order nets from existing higher 
' ^\ order nets and even classical (t,m, s)-nets. Such propagation rules for higher order 

P^ ■ nets were first considered by the authors in [9} and further developed in [2]. In 

^ . [1] Blokh and Zyablov established a very general propagation rule for linear codes. 

r^ ', This propagation rule has been extended to (t, m, s)-nets by Schiirer and Schmid 

C^ ' in [19]. In this paper we show that this propagation rule can also be extended to 

higher order nets. Examples indicate that this propagation rule yields new higher 

order nets with significantly higher quality. 
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^ '• 1 Introduction 

Quasi-Monte Carlo (QMC) rules are quadrature rules for approximating a high- dimensional 
rN ■ integral L -^^^ /(a;) da; by the average of certain function values, -^J2i=^ /(*«); where 

_c^_' Xq,Xi, . . . ,Xn_i are deterministically chosen integration nodes. The crucial question is 

how to choose the points Xi in order to obtain low integration errors. For this reason, 
several different classes of point sets and their properties in relation to integration prob- 
lems have been studied over the past decades (see, e.g., the monographs [HI [131 EO] for 
introductions to this topic). Note that by a point set we always mean a multiset, i.e., 
points are allowed to occur repeatedly. One prominent class of point sets commonly used 
in QMC algorithms are {t,m, s)-nets as introduced by Niederreiter (cf. [121 [131 [E] ) ■ It is 
known that (t, m, s)-nets can yield the optimal order of convergence (up to powers of the 
logarithm of the total number of points) of the integration error in QMC algorithms for 
functions of finite variation in the sense of Hardy and Krause (see again [T2l \T3\ E]). 
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In [5l |6], Dick generalised Niederreiter's (t,m, s)-nets to so-called higher order nets 
which have the convenient property that the corresponding QMC integration algorithms 
yield higher order convergence of the error for smoother functions. To be more precise, in 
[H] digital higher order nets were introduced, the construction of which, as it is also the 
case for Niederreiter's digital nets, is based on linear algebra over finite fields. In what 
follows, let ¥q be the finite field of prime-power order q. Moreover, we denote by N and 
No the set of positive integers and nonnegative integers, respectively. 

Definition 1. Letq be a prime power and letn, m, s G N. LetCi, . . . ,Cs benxm matrices 
over the finite field ¥q of order q. We construct g™' points in [0, 1)^ in the following way: 
For < h < q"^ let h = Hq + hiq + ■ ■ ■ + hm-iq"^^^ be the base q representation of h. 
Consider an arbitrary but fixed bijection rj : {0, 1, . . . , g — 1} -^ ¥q where 77(0) is the zero 
element in ¥g. Identify h with the vector h := (?7(/io), . . . , Tjijim-i)) G F"^- For I < j < s, 
we multiply the matrix Cj by h, 

Cyh7 =: (y,-i(/i),...,|/,-„(/i)) gF„ 

and set 

^h - ^ +■■■+ ^n 

Finally, set Xh := (xj^ , . . . , x)^ ) . The point set consisting of the points Xq, . . . , Xqm_i is 

called a digital net over ¥g. The matrices Ci, . . . ,Cs are called the generating matrices of 
the digital net. 

As it can be seen from Definition [H the properties of the points of a digital net (such 
as, e.g., their distribution in the unit cube) are determined by properties of the generating 
matrices Ci, . . . , Cg- These properties are, in the currently most general form of digital 
nets as introduced in [6], described by additional parameters t,a,(3, which is why those 
nets are referred to as (t, a, (3,n x m, s)-nets. The exact roles of the parameters t, a, and 
P are outlined in the following definition. 

Definition 2. Let n,m,a ^N, let < f3 < min(l, am/n) be a real number. Let ¥q be the 
finite field of prime power order q and let Ci, . . . ,Cs G F^^™" with Cj = {cj^i, . . . , Cj^n) ■ 
The digital net with generating matrices Ci, . . . ,Cs is called a digital {t, a, P,nx m, s)-net 
for a natural number t, < t < f3n, if the following condition is satisfied. For each choice 
of 1 < ij^y^ < • ■ ■ < ij^i < n, where Uj > for j = 1, . . . ,s, with 

il,l + \- il,mm{iyi,a} + h is,l + h is,mm{us,a} < /3n — t (l) 

the vectors 

are linearly independent overWg. 

Ift is the smallest non-negative integer such that the digital net generated by Ci, . . . ,Cs 
is a digital [t, a, /3, nxm, s)-net, then we call the digital net a strict digital {t, a, /3, nxm, s)- 
net. 

Note that Definition [2] implies that t must be chosen such that z/i -|- ■ ■ ■ + z/^ < m holds 
whenever ([T]) is satisfied. (Note that Uj < ij^i.) 



The definition of classical digital (t, m, s)-nets is obtained by choosing a = (3 = \ and 
m = n m Definition [2l 

Digital higher order nets are a subclass of general higher order nets which were intro- 
duced in [8]. While we denote digital higher order nets as (t, a, /3, ra x m, s)-nets, by which 
we emphasise the role of the generating matrices in the construction, general higher order 
nets are denoted as (t, a, /3, n, m, s)-nets. 

We give the definitions and some properties of (t, a, /3,n,m, s)-nets in base h. We 
follow [2] in our presentation. 

Let n, s > 1, 6 > 2 be integers. For i/ = (z/i, . . . , z/^) G {0, . . . , nY let \v\i = ^^^^^ i^j 
and define i^ = (ii^i, . . . , ii^ui, ■ ■ ■ , is,i, ■ ■ ■ , is,us) with integers 1 < ij^y. < ■ ■ ■ < ij^i < n in 
case Uj > and {ij^i, • • • , ij^uj } = in case Uj = 0, for j = 1, . . . ,s. For given u and v let 
a^ e {0,...,b- l}'"^'', which we write as a^ = (ai,i,_i, . . . , ai,j,_^^, . . . , a^,j^ i, . . . , a^^i^ ,^J. 

A generalised elementary interval in base 6 is a subset of [0, 1)'^ of the form 



s 


0— i 

u 




i=i 


"j,i=o 




«6{1,. 


.,n}\{ij,i,. 


.,b-,.,} 



From P[ Lemmas 1 and 2] it is known that for i^ G {0, . . . , nY and ij^ defined as above 
and fixed, the generalised elementary intervals J(i^, ajy) where a^ ranges over all elements 

from the set {0, . . . , 6 — l}'*^'^ form a partition of [0, l)** and the volume of J{iu, o,^) is 

6-Hi. 

We can now give the definition of (t, a, (3, n, m, s)-nets based on [SI Definition 4]. 

Definition 3. Let n, m, s, a > 1 be natural numbers, let < /3 < Tain{l, am/n) 
be a real number, and let < t < /3n be an integer. Let 6 > 2 be an integer and 
P = {xq, . . . , a;f,m_;^} be a multiset in [0, 1)''. We say that P is a (t, a, (3, n, m, s)-net in 
base b, if for all integers 1 < ij^,^. < ■ ■ ■ < ij^i, where < z/j < n, with 

s mm{uj ,a) 

j=i 1=1 

where for Uj = we set the empty sum Yli=i h,i ~ ^i ^^^ generalised elementary interval 
J(v, a^) contains exactly 6™^l^li points of P for each a^ G {0, . . . , 6 — 1}'^'\ 

A (t, a, /3,n, 771, s)-net in base b is called a strict {t,a, (3,n,m,s)-net in base b, if it is 
not a {u, a, (3, n, m, s)-net in base b with u <t. 

Informally we refer to (t, a, /3, n, m, s)-nets as higher order nets. 

Note that in the definition above the specific order of elements of a multiset is not 
important. The parameter t is often referred to as the quality parameter of the net. By 
the strength of the net one means the quantity f3n — t. 

The advantage of the more general concept due to [8] (in comparison to classical 
(t, -m, s)-nets) is that (t, a,/3, ?T,,m, s)-nets in base b can exploit the smoothness a of a 
function / (which is not the case for the classical concepts of (t, m, s)-nets and (t, s)- 
sequences). More precisely, we have the following result from [I] (and [SIE] for the digital 
case). 



Proposition 1. Let {xq, . . . , x^m-i} be a {t, a, /3, n, m, s)-net in base b. Let f : [0, 1]^ — )■ 
M have mixed partial derivatives up to order a > 2 in each variable which are square 
integrable. Then 



1 ^'""' 

"^[O'l]^ ^ h=0 



C)(5-(l-l/«)(/3n-«)(^^_i)«^) 



// {a;o, . . . , Xbm_i} is a digital {t, a, /3, n, m, s)-net in base b, then 



6™ — 1 



C(6-(^"-*)(/3n-t)"^). 



Additionally, the following results are known. If a = /3 = 1 and n = m, then the 
integration error is of order 0{b~"^~^^m^), see [T3] . 

Proposition [T] underlines the importance of knowing explicit constructions of higher 
order nets with a large value of (3n — t. 

In a series of papers, see for example [31 [151 HZl HI], so-called propagation rules for 
classical (t, m, s)-nets were introduced, which allow one to construct new digital nets 
from known ones and thereby improve on the parameters, in particular on the strength, 
of those nets. That such constructions are very useful can be seen in [12], where the 
best known parameters of classical (t, m, s)-nets are listed. A first step in establishing 
and using propagation rules for digital higher order nets (i.e., {t,a,(3,n x m, s)-nets) 
was made in the paper [9], where a series of propagation rules for such point sets were 
discussed. In the later paper [2], these propagation rules were extended to the case where 
the underlying nets need no more be digital, i.e., the rules apply to the more general class 
of {t, a, /3, n, m, s)-nets. 

In this paper, it is our aim to show a further propagation rule for higher order nets, 
which is a modification of a propagation rule for linear codes by Blokh and Zyablov 
[1]. Our main result, which will be prepared in Section [3] and stated in Section IH is an 
extension of a Blokh-Zyablov propagation rule for classical (t, m, s)-nets by Schiirer and 
Schmid in [TB]. Before (Section |5|), we are going to review duality theory of higher order 
nets, as it was introduced in [9] and [2]. Duality theory will be the main tool in proving 
our new propagation rule for higher order nets. Finally, we outline exemplary numerical 
results in Section |5l which are going to highlight the strength of our new propagation 
rule. 



2 Duality theory 

In this section we review the basics of duality theory for higher order nets, as this is going 
to be the key ingredient for showing our main results. Duality theory, as introduced for 
classical (t,m, s)-nets by Niederreiter and Pirsic [16], is a helpful tool in the analysis and 
construction of digital nets. In [£] it was extended to digital higher order nets, and it 
was further extended to the non-digital case in [2]. We keep close to [^ and [9] in the 
following two sections. 



2.1 Duality theory for digital higher order nets 

In [9], the dual of the row space of a digital (t, a, /3, n x m, s)-net was introduced. Given 
the generating matrices Ci, . . . ,Cs of a digital higher order net, let 



C=iCj\---\Cj)e¥: 
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The points of a digital net are obtained from the linear subspace C of F"* of dimension at 
most m, given as the row space of C by Definition [T] The dual space of C is given by 

AT = {x G F^^ : X ■ y = G F, for all y G P}, 

i.e., M is the null space of C . 

Let A = (ai, . . . , a^) G F^**, where a^ = (a^^i, . . . , aj^n) £ IF^. For 1 < j < s with 

aj 7^ let 1 < i/j < n denote the number of nonzero elements of a^ and let 1 < ij^^^ < 

■ ■ ■ < ij,i < n denote the indices of the nonzero elements aj^i^^^^, . . . , Oj-.j^.^ G F* = F^ \ {0} 

of aj (thus, Qj^i = for / G {1, . . . , n} \ {ij^i, . . . , ij,uj})- Let a > 1 be an integer. Then 

we define 

r for a,- = 0, 

I L/c=i h,k otherwise, 

and 

s 

i=i 

Let A/" C F"** and let \f/\ denote the number of elements of A/". For A, B G A/" we 
define the distance 

d„,„(A,B) = /x^,„(A-B). 

Further we define 

SaA'^) = min d„,„(A,B). 

A,BeAf 

We always have Sa,n{J^) > 1 and Sa,n{J^) > Sa',n{J^) for a > a' > 1. 

The following definition and result was first shown in [9j and is a generalisation of [16] . 

Proposition 2. The matrices Ci, . . . ,Cs & jpnxm g^ji^j-^f^ q digital (t, a, f3,n x m, s)-net 
over Fg if and only if 

5aA^) >Pn-t + l, 

where M is the dual space of the row space C. If Ci, . . . ,Cs generate a strict digital 
{to,a,P,n X m,s)-net overFq, where we assume that (in is an integer, then 

2.2 Duality theory for general higher order nets 

Here we review duality theory for higher order nets which also applies to point sets not 
obtained by the digital construction scheme. The basic tool are Walsh functions in integer 
base 6 > 2 whose definition and basic properties are recalled in the following. We repeat 
some background from [2] , where also the proofs of the results mentioned in this section 
can be found. 



Definition 4. Let 6 > 2 be an integer and represent A; G No in base b, k = Ha-ib"'^^ + 
■ ■ ■ + kq, with Kj G {0, . . . , 6 — 1}. Further let Ub = e'^'^^^'' be the 6th root of unity. Then 
the kth b-adic Walsh function walfc(x) : [0, 1) — t- {l,ujh, ■ ■ ■ ,1^^""*^} is given by 

walfc(x) = a;«^''"+-+«"^"-\ 

for X G [0, 1) with base b representation x = ^ib~^ + ^2&^^ H — ■ (unique in the sense that 
infinitely many of the ^i are different from 6 — 1). 

For dimension s > 2, x = (xi, . . . , Xs) G [0, 1)'', and k = {ki, . . . ,ks) G Ng, we define 
walfc: [0,l)^^{l,a;b,...,a;^i}by 

s 

walfc(a;) = ]^wa^.(xj). 
i=i 

The following notation will be used throughout the paper: By © we denote digitwise 
addition modulo 6, i.e., for x,y E [0,1) with base b expansions x = ^^i^z6"' and 
y — Xli^i Vib\ we define 

1=1 
where Q G {0, . . . , 6 — 1} is given by (i = C^i + rji (mod b). Let denote digitwise sub- 
traction modulo b (for short we use Qx := a;). In the same fashion we also define 
digitwise addition and digitwise subtraction of nonnegative integers based on the 6-adic 
expansion. For vectors in [0, 1)* or Nq, the operations and are carried out compo- 
nentwise. Throughout the paper, we always use the same base b for the operations and 
as is used for Walsh functions. Further, we call x G [0, 1) a 6-adic rational if it can be 
written in a finite base b expansion. The following simple properties of Walsh functions 
will be used several times. 

For all A;, / G No and all x,y E [0, 1), with the restriction that if x,y are not 6-adic 
rationals, then x y is not allowed to be a 6-adic rational, we have walfc(x) • wali{x) = 
walfcez(a;) and walA;(a;) ■ walfc(y) = walfc(x y). Furthermore, walfc(a;) = waleA:(a;). 

Now we turn to duality theory for nets. Let K-l^ = {0, . . . ,b^ — ly. We also assume 
there is an ordering of the elements in /C^ ^ which can be arbitrary but needs to be the 
same in each instance, i.e., let JCl^. = {k^, . . . , kbsr_i}. (Note that j/C^^I = b^^.) By this 
we mean that in expressions like J^keK" ' i^k,i)k,ieic'' ) and {ck)keic= the elements k and 
/ run through the set /C^ ^ always in the same order. 

The following b^^ x 6**" matrix plays a central role in the duality theory for higher order 
nets, 

W.:=(wal,(6-Z)),^^^^.^. 

We call Wr a Walsh matrix. 

In the following we denote by A* the conjugate transpose of a matrix A over the 
complex numbers C, i.e.. A* = A . The following lemma was shown in p]. 

Lemma 1. The Walsh matrix^ ^ is invertible and its inverse is given byyV~^ = b^'^'^W*. 

Let now 6 > 2 and r, A^ > 1 be integers. For a multiset P = {xq, . . . , x^-i} in [0, 1)* 
and k E K-l^ we define 

N-l 

Ck = Ck{P) := y^walfc(a;fe) 

h=0 



(note that |cfc| < A^ and cq = N) and the vector 

(7 = (7(P):=(c,W^.^. (2) 

For a = (ai, . . . , a^) G /C^^ define the elementary 6-adic interval 

We also have the following lemma. 
Lemma 2. We have 



^ walfc(a; Q y) = I 



|/C^_J ifx,yeEaforsomeae}Cl,., 
otherwise. 



Let X G Ea for some a G /C^^. Then, using Lemma [21 we have 

Tjt:—^ ^ Cfcwalfc(a;) = — — ^ y^walfc(a;fe Q x) 

iV-l ^ 

= y^ jiTrT y^ walfc(a;/i ea;) 

= \{h : x/i e -Ea}| =: rria. 

Definition 5. Let b > 2 and r, A^ > 1 be integers. Let P = {xq, . . . , xi\i_i} be a multiset 
in [0, 1)" and let JCl^ = {0, . . . , 6^ - 1}^ 

1. For a e }Cl^ let 

rria = ma{P) := \{h : Xh G -Ea}| 

and 

M = M(P) := (m„)„,^.^ . 

Then we call the vector M the point set vector (with resolution r). 

2. The vector C = C{P) from ([2]) is called the dual vector (with respect to the Walsh 
matrix W^). 

3. The set 

Vr = Vr{P):={kelCl,:Ck^Q} 

is called the dual set (with respect to the Walsh matrix W^). 

The relationship between a point set vector and its dual vector is stated in the following 
theorem. 

Theorem 1. Let P = {xq, . . . ,a;7v-i} be a multiset in [0, 1)'^ and let r G N. Let M be 
the point set vector with resolution r and C be the dual vector with respect to W^ defined 
as above. Then 

^ WrC = M and C = W;M. (3) 



I A"* 



The vector C carries sufficient information to construct a point set in the following 
way: Given C, we can use Theorem [1] to determine how many points are to be placed 
in the interval Ea, a G /C^r- ^^ remark that at the functional level, the vector C is 
comparable to the generator matrices of a digital net, which completely determine the 
point set. 

Note that for the {t,a, f3,n,m,s)-net property it is of no importance where exactly 
within an interval Ea, a G /C^„, the points are placed. Hence we can reconstruct a net 
from a dual vector with respect to W^ provided that r > \_f3n\ — t. In words, if one knows 
the dual vector of a net, then one can use this dual vector to obtain the net via Theorem [T] 
provided that the resolution is greater than or equal to the strength of the net. 

In analogy, the dual space of a digital net also allows us to reconstruct the original 
point set, see [16|. Although C is different from the dual space for digital nets, it contains 
the same information and can be used in a manner similar to the dual space. In case P 
is a digital {t,a,(3,n x m, s)-net, the dual set Vn defined in Definition [5] coincides with 
the dual space defined in [S] intersected with /C^„, and if P is a digital (t,m, s)-net, it 
coincides with the dual space in p]5] intersected with /C^^. 

Although the above results hold for arbitrary point sets, in the following we consider 
point sets which are nets and show how to relate the quality of a (t, a, /3, n, m, s)-net to its 
dual set. To this end we need to introduce a function which was ffist introduced in [6] in 
the context of applying digital nets to quasi-Monte Carlo integration of smooth functions 
and which is related to the quality of suitable digital nets. For fc G No and a > 1 let 

(i.\ - j "1 "^ ^ «min(i/.a) for A; > 0, 

^"^'^^ ~ \ for fc = 0, 

where for fc > we assume that k = Kib""^'^ + ■ ■ ■ + Kj.b""'"^ with < ki, . . . , K,,y < b 
and 1 < a^ < ■ ■ ■ < ai. Note that for a = 1 we obtain the well-known Niederreiter- 
Rosenbloom-Tsfasman (NRT) weight (see, for example, [IT, Section 7.1]). 

For a vector k = {ki, . . . ,ks) G Nq we define HaiM) = Haiki) + ■ ■ ■ + fia{ks) and for a 
subset Q of ICl^ with Q \ {0} ^ and a > 1 define 

Pa(Q) := , min Ua(fc). 

keQ\{0} 

For Q C {0} we set p„(Q) = r + 1. 

Let P = {xq, . . . ,xn-i} C [0, 1)''. In the following we consider for which cases we 
have Vr{P) = {0} (note that G Vr{P) for any point set P with at least one point). If 
Vri^P) = {0}, then we have Cq 7^ and Cfc = for all fc G /C^^ \ {0}. By Theorem [1] we 

have M{P) = cob~'''^{l, 1, . . . , 1)^, that is, each box Ea contains exactly cob^^"^ points for 
all a G /C^ ^ and P consists of A^ = cq points altogether. This is the only case for which 

VriP) = {0}. 

Conversely, since the number of points in Ea must be an integer, it follows that 
Cob^^"^ G N, i.e., b^'^ divides Cq and therefore 6^* divides A^. From this we conclude that 
if we choose a resolution r G N such that 6'"'^ > A^, i.e., r > -log^ A^, then Vri^P) 7^ {0}. 
For a net with N = b"^ points this means that we require r > m/s. 

The following theorem establishes a relationship between Pa{Q) and the quality of a 
(t, a, (3,n,m, s)-net. 



Theorem 2. Let P = {xq, . . . , Xf)rn_i} c [0, 1)'' be a multiset. Then P is a {t, a, /3, n, m, s)- 
net in base b if and only if Pa(J^ii3n\-t) > [f3n\ —t + 1. If P is a strict (to, a, (3, n,m, s)-net 
m base b, then Pa(Pipn\~to) = L/^^J - to + 1- 

Let now P be a strict (to, a, /3, ra,m, s)-net in base b. Let r > \_Pn\ — to- Then 

'Dr ^ ^L/3"J-to and P, \ ©L/^nJ-to ^ ^,r \ ^llM-to- ^^^ ^^^ ^ ^ ''^b-r- \ ''^U/JnJ-to ^^ ^^^^ 

Pa{k) > [(3n\ — to + 1. Theorem [2] imphes that Pa(^[;3nj-to) = L/^'^J ~ ^o + 1 and hence 
PaiJ^r) = Pa{T^[i3n\~to) = [/^^J — to + 1- In particular, for aU r,r' > [f3n\ — to we have 

paiVr) = PaiVr') = p„(2)„) = [M - to + 1, (4) 

since n > [f3n\ — to. 

3 A generalized Blokh-Zyablov construction 

In this section we generalise the Blokh-Zyablov construction for codes from [4j and for 
classical nets from [18]. The proofs closely follow those in [18] . 

We call A/" C F^"" an {{s,n),a,K,6)g-spa:Ceii K = \Af\ and 6 = 6a,n{M). If A/" is a linear 
subspace of F^* then we call A^ a linear [s, n, a, fc, 5]q-space if jA/"! = q^ and 5 = Sa,n{-^)- 

Note that an ((s,n), 1, ii', (5)g-space is a generalized {{s,n), K, 6) g-code and a linear 
[(s, n), 1, k, 6]q-spsice is a generalized linear [(s, n), k, 6]q-code in the sense of Schiirer and 
Schmid [H], [19]. Note, furthermore, that in the special case n = 1 an [(s, 1), k, 6]q-code is 
an ordinary linear [s , k , 6]q-code (see [18] or [19] for details). 

Lemma 3. Let a,a' > 1 and let 

{0}=AfoCAfic---C M' = F; ■ ■ 



Q 



be a chain of linear [s' , n' , a' , k'^, 6'^]q-spaces and let \i, . . . ,\n's' € F^ '^ denote vectors 
such that Af^ is generated by yi, . . . , v^^ for u = 1, . . . , r. We call the spaces N^ the inner 
spaces. 

Let Mu denote (not necessarily linear) {{s,n),a,Ku,Su)q^-spaces with g„ = g*^" and 
Cu = k'^ — k'^_^ for u = 1, . . . ,r. We call the spaces Mu the outer spaces. 

Let M be given by 

A/" = < ^ 0«(x„) : x„ G A/'„ /or u = 1, . . . , r l , (5) 

where (pu '■ IF"* — )■ F^'" ** replaces each symbol (regarded as a vector of length Cu over 
¥q) of an element from Mu by the corresponding linear combination of the e„ vectors 
Vfc' _ +1, . . . , Vfc/^. The elements in the resulting vector are grouped such that column (a', r') 
(with 1 < a' < s' and 1 < r' < n' ) from M^ and column (a, r) (with 1 < a < s and 
1 < T < n) from Mu determine column {{a — l)s' + a', (r — l)n' + r') in the resulting 
element in M. 
Then M is an 

{{ss',nn'),aa',Ki- ■ ■ Kr, min 5„5^)„ — space. 

l<«<r,|A/'„|>l 



Proof. The proof follows along the same lines as the proof of [19^ Theorem 3]. We have 
A/" C Fg" '^^ by the definition of the mappings 0„. The number of elements in JV is given 
hj Ki- ■ ■ Kr which follows from ([5]) and the fact that all elements given by X]n=i 0m (^n) 
for Xu G Mu and w = 1, . . . , r are distinct, which we show in the following. 

Let X, y be two elements M and let Xij, y„ G ^4, 1 < m < r, denote the elements 
defining x and y using the sum in ([5]). If x„ = y^ for 1 < m < r, then it follows that 
X = y. Assume now that x 7^ y, then there exists a largest integer 1 < m* < r such that 
Xm* 7^ Yu*- This implies that |A/'t(*| > 1. 

Since 0„ is linear, we have 

u* — l r 

x-y = ^ (/)«(x„ -y„) + (/)^,*(x„. -y„.) + ^ </)n(x„-y„), 

«=1 u=u*+l 

where the first sum is a concatenation of elements from Ml*_i and the second sum is 
(since x^, = y^ and the linearity of 0„ for u* < u < r). Further, (/)„*(xu* — y„*) is a 
concatenation of elements from M^* \J^u*-i- ^^ denote the elements in this concatenation 
by w-^''^ for 1 < j < s and 1 < k < n. 

Thus, X — y is a concatenation of elements from A/'^. \ A/^*_i, since the sum of an 
element from Af^*-! ^^^ ^^ element from M^^, \M^*_i is again an element of A/^. \ A/'^.-i- 

Let 

j=l k=l 

where Uj denotes the number of nonzero elements and Zj ^ are the positions of the nonzero 
elements in x„. — y^.. Furthermore, let 





s' ™"(°''^kj.fc) 


fia'^n'i^^''"^) 


= Z^ Z^ '^'j',k',j,kJ 




j'=l k'=l 



where i^^j^k denotes the number of nonzero elements and i'j>^kij^k denotes the positions of 
the nonzero elements in w^-^''^''. Then we have 

Waa',nn'(,X, y ) = Haa' ,nn'\^ Y) 

s s' min(a,iyj)mm{a',;y^,, ^. J 

> Y" \^ \^ V " I ^^^'^ - 1)^ + V,fc'j,fc if ^i' ^r > 

~ -^ -^ -^ -^ 1 otherwise 

j=i j'=i fe=i fc'=i >- 

j=l j'=l fc=l fc'=l 

= / . / . '^hk 2_^ 2^ '^'j',k',j,k 

j=l k=l j'=l k'=l 

>5a,ni^u*)Sa',n'iK*)- 

The result thus follows. D 
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Lemma 4. Let {0} = Mq d M[ d ■ ■ ■ d M'^ = F^' he as in Lemma\M with a' = n' = 1. 
Let Vi, . . . , Vs' G Fg be defined as in Lemma\^ where we assume without loss of generality 
that the first i — 1 elements of Vj are 0. This means that (vi, . . . , Vg/)^ forms an s' x s' 
upper triangular matrix. 

Moreover, let qu and Cu be defined as in Lemma 0. For positive integers si < S2 < 
■ ■ ■ < Sr let Nu, M = 1, . . . , r, be a (not necessarily linear) ((s^, n), a, Ku, 5u)q^-space. 

Then we can construct an 

((siCi + ■ ■ ■ + SrCr, n),a,Ki- ■ ■ Kr, mill SuS'^^q — space. 

l<M<r 

Proof. The construction is analogous to the construction in the proof of [19^, Theorem 4] 
and follows in three steps. 

• Let s = Sr and construct new spaces Aiu by embedding each space A4 in ^qu ■ 
This is achieved by prepending G Fg^"'^"'"'^ to each element of A/^. 

• Use Lemma [3] with A/'g, . . . , A/"/ as inner spaces and A^i, . . . , J^r as outer spaces. 
This yields a new space A^. 

• Finally, for i = 1, . . . ,s', choose u minimal such that the condition ei + ■ ■ ■ + e„ > i 
holds and construct a space Af from A/f. This can be done by deleting all columns 
[k, t) in M with k = Os' + i,ls' + i, . . . ,{s — Su — l)s' + i. 

The total number of deleted blocks is ei(s — Si) + ■ ■ ■ + er{s — Sr), so the length of 
M is SiCi + ■ ■ • + SrCr- The deleted positions are G F^ for each 0„(x„), either due 
to a in Vj or due to a appended to x„ G Mu- Note that this procedure neither 
influences the dimension nor the weight of A4 . 

D 



4 A Blokh-Zyablov propagation rule 

4.1 A Blokh-Zyablov propagation rule for digital higher order 
nets 

Lemma m can be applied to digital higher order nets which yields a new propagation rule. 
This propagation rule generalises the Matrix-product construction in [9]. As there were 
Propagation Rules I-XIV in [9], we call our new rule Propagation Rule XV. 

Theorem 3 (Propagation Rule XV). Let vi, . . . , v^/ and M^, qu and Cu for u = 1, . . . , r 
be defined as in Lemma\^ Let a > 1. For positive integers Si < S2 < ■ ■ ■ < Sr let Vu, 
1 < u < r, denote digital (t„, a, f3u, Uu x m„, Su)-nets over F^^. 

Then a digital {t,a,P,n x m,s)-net overWg can be constructed, where 



m = 

u=l 



r 
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E' 



n = y e,,n. 



■U"'U1 



M=l 



/3 = min(l, am/n), 
t < (3n + l- mill (/3„n„ - t„ + 1)5^. 

l<M<r 

Proof. Let A/'u denote the dual space of Vu for u = l,...,r. Then A/'^ is a hnear 
[(s,n), Q;,m„, /3„n„ — t„ + l]q^-space. The result then follows by applying Lemma H] us- 
ing A/'o, . . . , A/"/ as inner spaces and A/i, . . . ,Mr as outer spaces to obtain a space M and 
applying Proposition [2J D 

4.2 A Blokh-Zyablov propagation rule for higher order nets 

In this subsection we consider higher order nets which are not necessarily digital. Again 
we have the same theorem as in the previous case by considering the dual set instead 
of the dual net. Applying Lemma S] and Theorem [2] yields the following result, which is 
Propagation Rule 17 for general higher order nets (as there were Propagation Rules 1-16 
for general higher order nets in [2]). 

Theorem 4 (Propagation Rule 17). Let vi, . . . , vy and M!^, g„ and e„ for u = 1, . . . , r he 
defined as in Lemma^ Let a > 1. For positive integers Si < S2 < ■ ■ ■ < Sr let Vu denote 
{tu, (3u, a, Uu X m„, Su)-nets in base hu = fe"^" for 1 <u <r. 

Then a {t,a,l3,n x m,s)-net in base b can be constructed, where 

r 
u=l 



m = 

u=l 



u=\ 

T 



n = 

u=l 

/3 = min(l, am/n), 
t < /3n + l- min (/3„n„ - t„ + 1)6'^. 

l<u<r 

Remark 1. Note that in Lemma\^ Lemma^ Theorem\^and Theorem^it is sufficient 
to require A/"/ C F^ instead of M'. = F^ . Indeed, this slight modification is possible as one 
can use a trivial space Mu = {0} as an outer space in Lemma\^ which does not influence 
the minimal distance of the newly obtained space (cf. Remark 4 in fTEP - 

5 Numerical results 

Theorems [3] and H] provide previously unknown propagation rules for higher order nets. 
To show how powerful these are, we state exemplary numerical results for the digital 
case, but it should be noted that the results are also true for the more general case, as the 
propagation rules yield the same parameters and we always start from existing digital nets, 
which are by definition a subclass of general higher order nets. In Tables [TH31 we present, 
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for selected values of m and s, the results obtained when we use different propagation 
rules for the case q = 5, a = 2, and /3 = 1. As we restrict ourselves to considering only 
digital nets, we only refer to propagation rules with roman numbers within this section. 

Table [T] covers the case where s = 5, Table [2] the case s = 15 and Table [3] the case 
s = 25. In all tables we consider m between 15 and 30. Since a and (3 are fixed, and 
different propagation rules might yield different ratios of n and m, it is most useful to 
compare the strengths of the nets obtained. As outhned above, the strength of a digital net 
refers to the value of a = cr(/3, n, t) = [in — t. Note that, in some of our new propagation 
rules, one can make many different choices of smaller nets that might yield a bigger net 
with the same parameters. We only give the best values of the strength a we can obtain 
by going through a number of possible choices of the smaller nets involved. 

Our tables below can be seen as an extension of Table 3 in [9|, Section 4]. To be more 
precise, we consider the following quantities in Tables [1] and [2l 

• o"dir: The strength of a digital (t, 2, 1, 2m x m, s)-net over F5 using the generating 
matrices of an existing classical digital (t', m^ 2s)-net over F5, where we then obtain 

(cf. m 

t < 2min|m,t'+ - | . (6) 

We call this construction method the direct construction method. 

• (Tvii: The strength of a digital net constructed from a digital (ti, 2, 1, 2mi x mi, si)- 
net Pi and a digital (t2,2,l,2m2 x m2,S2)-net P2 over F5 using Propagation Rule 

VII in [0], where Pi and P2 are obtained by the direct construction method from 
classical nets. Here, n = 2m. 

• o"viii: The strength of a digital net constructed from a digital (ti, 2, 1, 2mi x mi, si)- 
net Pi and a digital (^2, 2, 1, 2m2 x m2, S2)-net P2 (si < S2) over F5 using Propagation 
Rule VIII in ^, where Pi and P2 are obtained by the direct construction method. 
Here, n = 2m. 

• aix'. The strength of a digital net constructed from a digital (ti, 2, 1, 2mi x mi, si)- 
net Pi, a digital (^2, 2, 1, 2m2xm2, ■S2)-net P2, and a digital (ta, 2, 1, 2m3xm3, S3)-net 
-P3 (si < S2 < S3) over F5 using Propagation Rule IX in P, where Pi, P2, and P3 
are obtained by the direct construction method. Again, n = 2m. 

• axi- The strength of a digital net obtained by using Propagation Rule XI from 
[9] with r = 2, where the higher order nets plugged into Rule XI are obtained 
by the direct construction method from classical nets. Note that this rule can, 
since r = 2, only be applied for the cases where m is even, and we then have 
n = m,. As Propagation Rule XI with r = 2 only yields higher order nets with 
an even value of s, we obtain the values in the Tables [TH3] by projection from 6- 
dimensional, 16-dimensional, and 2 6- dimensional nets, respectively (this is allowed 
due to Propagation Rule V in |9]). 

• cTxi+viii: The strength of a digital net obtained by first applying the direct con- 
struction method to classical nets, then applying Propagation Rule XI with r = 2 
in |9], and then using pairs of the newly obtained nets to apply Propagation Rule 

VIII in mj . Note again that we need to restrict ourselves to even cases of m, here, 
as we first apply Rule XI with r = 2. For o"xi+viii, we again have n = m. 
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m 


Cdir 


CTVII 


fviii 


ClX 


Cxi 


CXI+VIII 


cxv 


15 


24 


12 


17 


14 








16 


26 


14 


18 


16 


14 


9 


19 


17 


28 


14 


20 


17 








18 


30 


16 


21 


18 


16 


10 


21 


19 


32 


16 


22 


20 








20 


34 


18 


24 


20 


18 


12 


25 


21 


36 


18 


25 


21 








22 


38 


20 


26 


22 


20 


13 


27 


23 


40 


20 


28 


24 








24 


42 


22 


29 


25 


22 


14 


29 


25 


44 


22 


30 


26 








26 


46 


24 


32 


26 


24 


16 


33 


27 


48 


24 


33 


28 








28 


50 


26 


34 


29 


26 


17 


35 


29 


52 


26 


36 


30 








30 


54 


28 


37 


32 


28 


18 


37 



Table 1: a-values depending on m (15 < m < 30) for a = 2, /3 = 1, g = 5, and s = 5. 

• axv'- The strength of a digital net by first applying the direct construction method 
to classical nets, and then Propagation Rule XV with r = 2. (Theorem [3]). As the 
inner spaces (see the proof of Theorem [31 and Lemma H]) we use (extended) Reed 
Solomon Codes over F5. These codes are linear [5,j, 6 — j]-codes over F5 (see |19j). 
In view of Remark [T|, we can use Reed Solomon Codes with parameters [5,0,6]5, 
[5, 2, 4)5, and [5,4, 2)5 as the inner codes, hence r indeed equals 2 and ei = 62 = 2 
in Theorem |3l which is why we again need to restrict ourselves to even cases of m 
here (in general one can of course also obtain odd values for m by different choices 
of Cj and rrii). In the language of Theorem [31 we have 61 = 4 and ^2 = 2 in this case. 
For o"xv we have n = 2m, and, as for cxxi, we obtain the values in the Tables [TH3] by 
projection from 6-, 16-, and 26-dimensional nets, respectively (again by Propagation 
Rule V in |9]). 

The results in Tables [TH3] show that our new Propagation Rule XV is superior to the other 
propagation rules tested for dimensions 15 and 25, and that also the combination of other 
rules, such as that of Rules XI and VIII, does not yield better results. In dimension 5, 
though still competitive. Propagation Rule XV is outperformed by the direct construction 
method. It is likely that the direct construction performs better for lower dimensions as 
the t-values of classical digital nets are very small in this case. Furthermore, the bound 
(|6]) depends critically on the value of s, so it is natural that lower dimensional results 
for the direct construction method are stronger. The fact that the direct construction 
method works well for low dimensions is in line with numerical results in [10], where a 
similar phenomenon was observed. 

We emphasise that our examples are just illustrations and can by no means system- 
atically cover all cases one might theoretically consider; to be more precise, we have the 
following restrictions in Tables [TH3l 

• We only show particular choices of the parameters involved. We restrict ourselves 
to some illustrative cases. 
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m 


""dir 


o'vii 


"■viii 


o^ix 


o'xi 


fXI+VIII 


o'xv 


15 





4 


5 


6 








16 


2 


4 


5 


8 


8 


6 


13 


17 


2 


4 


6 


8 








18 


4 


6 


8 


9 


10 


8 


17 


19 


4 


6 


9 


10 








20 


6 


6 


10 


12 


12 


9 


19 


21 


8 


8 


12 


13 








22 


8 


8 


13 


14 


14 


9 


21 


23 


10 


10 


13 


14 








24 


12 


12 


14 


16 


16 


12 


25 


25 


14 


14 


16 


17 








26 


16 


16 


17 


18 


18 


13 


27 


27 


18 


18 


18 


20 








28 


20 


20 


20 


20 


20 


14 


29 


29 


22 


22 


22 


22 








30 


24 


24 


24 


24 


22 


16 


33 



Table 2: o"- values depending on m (15 < m < 30) for a; = 2, /3 = l,g = 5 and s = 15. 



m 


""dir 


CTVII 


CTVIII 


o^ix 


c^xi 


c^xi+viii 


o^xv 


15 










2 








16 










2 


4 


2 


5 


17 










2 








18 










2 


6 


4 


9 


19 










2 








20 










2 


8 


5 


11 


21 










4 








22 










4 


10 


6 


13 


23 





2 


2 


5 








24 





2 


2 


5 


12 


8 


17 


25 





2 


4 


6 








26 





2 


4 


6 


14 


9 


19 


27 


2 


4 


5 


8 








28 


2 


4 


5 


8 


16 


10 


21 


29 


4 


4 


5 


9 








30 


6 


6 


6 


10 


18 


12 


25 



Table 3: a-values depending on m (15 < m < 30) for a = 2, P = 1, q = 5, and s = 25. 
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• We do not consider all possible combinations of different propagation rules, as this 
would lead to a too high number of parameters. 

• Not all propagation rules are applicable for all sets of parameters. This is indi- 
cated by void cells in the tables in cases where a certain propagation rule was not 
applicable. 
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